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=  voltage  exchange  parameter  =  V^T/  rh  )2 l(V ,y¥/Jlj  =  ®(l-p)% 

=  mass  exchange  parameter  =  ( m  =  (l-r)/% 

=  exit  mass  fractions  of  Xe°,Xe+1,Xe+2,Xe+3,  where  f0  +  f,  +  f2  +  f3  =  1 
=  ion  mass  fraction  at  exit  =  f]+f2+f3 

=  reduced  ion  mass  fractions  at  exit,  f  i+f  2+f  3=1.  where  f  i=fi/fi,  f  2=f2/f,,  f  3=f3/fi 
=  Faraday  constant,  96,485  coulombs/mol  of  charge 
=  velocity  distribution  function  of  ions  and  neutrals 
=  Earth’s  gravitational  constant  at  sea  level,  9.806  m/s2 
=  anode  discharge  current 
=  integrated  beam  current  =  f,  Q  rii  (&/JI) 

=  integrated  axial  beam  current 

=  specific  impulse,  independent  of  unit  system,  Isp  =  T  /  m 
=  charge  density  in  the  plume  measured  from  a  Faraday  probe,  A/sr 
=  anode  propellant  mass  flow  rate 

=  molecular  weight  of  propellant  (xenon  =  0.1313  kg/mol) 

=  input  power  to  anode  =  Va  Ia 
=  jet  power 

=  average  charge  of  propellant  ions  in  a  trimodal  jet=  f  ,  +  2f  2  +  3f  3 
=  unit  of  charge,  1.609  x  10"19  Coulombs 
=  fractional  anode  electron  current,  electron  recycle  fraction 
=  fractional  ion  beam  current,  current  utilization  efficiency  =  Ifieum  /  Ia 
=  minimum  fraction  of  electron  current  required  to  sustain  discharge 
=  thrust 

=  applied  anode  potential 
=  exit  velocities  of  Xe°,Xe+1,Xe+2,Xe+3 
=  exit  velocity  of  ions  and  neutrals 
=  exit  velocity  of  ions 

=  reduced  velocity  ratios,  where  y0  =  v(/l \'i I,  y,  =  v  ,/l V| I,  y2  =  v2/lvil,  y3  =  vVI  V|  I 
=  indicator  of  the  ‘quality’  of  ionization  in  ®vdf.  «  =  cl)vDi  /fi  ~  1.10-0.10Q 
=  relationship  between  particle  momentum  and  current  at  azimuthal  position  0,  y  ~  («/Q)1/2 
=  fractional  loss  of  acceleration  voltage 
=  voltage  utilization  efficiency  =  AV  /  Va 

=  average  ion  acceleration  voltage  =  (f  i  5Vi+  2f  2  §V2+  3f  3  5V3)  /  Q 
=  acceleration  potential  of  Xe+1,  Xe+2,  Xe+3 

=  ionization  potential  of  Xe+1,  Xe+2,  Xe+3  (12,  33,  65  eV  from  neutral  ground  state) 

=  energy  efficiency  =  Pjet/Pin 
=  experimental  thruster  efficiency 

=  azimuthal  position  from  thruster  centerline,  0=0°  on  thrust  axis 
=  plume  momentum  divergence  half-angle,  X=0°  on  thrust  axis 
=  propellant  utilization  efficiency,  ®  =  ®vdf  ^div 

=  velocity  distribution  loss  component  of  propellant  utilization,  4>vdf  =  «  fi 
=  divergence  loss  component  of  propellant  utilization,  ®div  =<cos(0  )>,m2  =  <cos(0)>/2 
=  output  moles  of  charge  per  input  moles  of  propellant,  average  particle  charge  =  f;  Q 


<  >m  =  mass  weighted  average  quantity  in  the  plume  (-7i/2<0<7t/2) 

<  >„,v  =  momentum  weighted  average  quantity  in  the  plume  (-7t/2<0<7t/2) 

<  >j  =  charge  flux  weighted  average  quantity  in  the  plume  (-7t/2<0<7t/2) 

(0)  =  variable  evaluated  at  azumuthal  position  0,  units  of  sr"1 

z  =  ion  charge  state  =  1,  2,  3  for  Xe+1,  Xe+2,  Xe+3 
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Introduction 


Standardization  of  experimental  methods,  facilities,  and  apparatus  in  the  electric  propulsion  community  has 
previously  been  proposed. 1  A  standard  approach  to  Hall  thruster  efficiency  analysis  is  also  necessary  to  enable 
consistent  comparisons  of  thruster  performance.  Present  Hall  thruster  performance  models  arbitrarily  set  the  anode 
efficiency  equal  to  the  product  of  several  utilization  efficiencies,  including  combinations  of  current,  voltage,  charge, 
mass,  and  propellant  utilization.2'3,4'5  This  method  will  often  successfully  interpret  experimental  measurements,  but 
it  is  not  based  on  first  principles  and  there  is  no  assurance  that  it  accurately  interprets  global  plasma  properties  or 
relevant  loss  mechanisms. 

A  standard,  logical  architecture  is  necessary  for  the  analysis  of  Hall  thruster  performance  data,  telemetry 
(Va,  Ia,  and  m  ),  and  plume  measurements  that  is  based  on  a  consistent  set  of  definitions,  energy  conservation,  and 
Newton’s  Second  Law.  In  electrostatic  propulsion,  thrust  is  the  result  of  the  utilization  of  propellant,  applied 
acceleration  potential,  and  discharge  current.  It  is  therefore  reasonable  to  express  thrust  efficiency  as  a  product  of 
the  utilization  of  these  three  elements. 

In  this  paper,  a  Hall  thruster  performance  model  is  derived  by  analytical  separation  of  thrust  efficiency  into 
the  product  of  energy  efficiency  and  propellant  utilization  efficiency,  which  are  less  than  unity  under  all  operating 
conditions.  This  simple  framework  reconciles  various  usages  of  propellant  utilization  and  exposes  the  distinction 
between  energy  efficiency  and  thrust  efficiency.  Energy  efficiency  is  expressed  as  a  product  of  voltage  and  current 
utilization  efficiencies,  and  is  rigorously  separated  from  losses  associated  with  the  jet  vector  properties.  The  model 
accounts  for  the  effects  of  the  ion  velocity  distribution  function  (VDF),  ionization  fraction,  multiply  charged  ion 
species,  electron  recycle  fraction,  momentum  plume  divergence,  and  the  loss  of  acceleration  potential  within  the 
discharge  chamber. 

The  proposed  efficiency  architecture  is  compared  to  past  methodologies,  with  a  historical  perspective  of 
early  ion  engine  analyses  and  similar  efforts  from  the  former  Soviet  Union.  To  illustrate  the  fidelity  of  the  proposed 
analytical  model  and  to  demonstrate  the  utility  of  various  diagnostics  in  determining  Hall  thruster  utilization 
efficiencies,  experimental  performance  and  plume  measurements  from  a  single  operating  condition  of  a  laboratory 
Hall  thruster  is  presented.  The  nature  of  analytical  factorization  allows  for  the  characterization  of  loss  mechanisms 
and  interpretation  of  causal  plasma  phenomena,  but  does  not  account  for  parasitic  losses  such  as  cathode  flow  and 
magnet  power. 
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Historical  Perspective  and  Recent  Efforts 


Investigations  of  ion  acceleration  and  electron  transport  using  Hall  thruster  technology  began  in  the  early 
1960s  in  the  United  States  and  former  Soviet  Union.6,7’8’9’10’11'12'13’14’15'16  The  focus  of  US  electric  propulsion 
research  shifted  primarily  to  ion  engine  technology  in  the  early  1970s17,  whereas  investigations  in  the  USSR 

18  19  20 

continued  Hall  thruster  advancements  throughout  the  following  decades.  '  ’ 

Analytical  factorization  of  ion  engine  thruster  efficiency  was  described  as  early  as  1975  by  Masek  et  al.  in 
a  review  of  ion  engine  performance.21  Thrust  efficiency  was  factored  into  the  product  of  energy  efficiency  and 
propellant  utilization  efficiency  using  Newton’s  Second  Law,  with  terms  accounting  for  losses  due  to  multiply 
charged  ions  and  beam  divergence.  However,  the  derivation  was  condensed  and  justification  for  treatment  of  beam 
divergence  and  multiply  charged  ions  was  not  presented.  The  methodology  was  not  adopted  or  cited  in  the  ion 
engine  community  and  is  absent  in  modern  analysis. 

In  the  1990s,  the  manifestation  of  Russian  Hall  thruster  technology  in  the  Western  spacecraft  community 
catalyzed  a  resurgence  of  Hall  thruster  research.  '  ’  ’  Technology  transfer  that  began  in  the  early  1990s 
brought  invaluable  benefits  and  advancements  to  Western  Hall  thruster  development,  but  much  of  the  earlier  Soviet 
progress  in  Hall  thruster  research  that  was  published  in  the  Russian  language  was  not  translated.  The  limited 
availability  of  translated  documents  from  this  extensive  literature  inevitably  impeded  the  transmission  of  knowledge 
and  progress.  As  a  result,  several  important  contributions  from  Russian  research  have  not  been  widely  disseminated 
in  the  West,  including  the  analytical  factorization  of  thruster  efficiency. 

Factorization  of  Hall  thruster  efficiency  was  outlined  in  a  manual  on  stationary  plasma  engines  by  Belan, 
Kim,  Oranskiy,  and  Tikhononv  from  the  Kharkov  Aviation  Institute  in  1989. 28  This  seminal  document  was  later 
referenced  and  the  performance  methodology  summarized  by  Bugrova  et  al.  without  explicitly  stating  the  nature  of 
the  earlier  derivation.29  Kim’s  highly  cited  1998  paper  on  processes  that  determine  Hall  thruster  efficiency30  and 
other  contemporary  publications'  '  '  continued  the  elaboration  of  the  Kharkov  Aviation  Institute  methodology. 
These  approaches  showed  similarities  to  the  analysis  presented  here,  and  demonstrated  relationships  between 
experimental  variables  and  plasma  phenomena  that  affect  thruster  performance.  Over  time,  factorized  performance 
analyses  were  substituted  with  thruster  anode  efficiency  as  the  product  of  appointed  utilization  efficiencies  that 
capture  Hall  thruster  losses  as  a  function  of  various  physical  processes  with  varying  levels  of  completeness.  '  ’  '  ’  "  ’ 

35,36,37,38 
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A  textbook  by  Grishin  and  Leskov  analyzed  thruster  performance  based  on  energy  efficiency  and  thrust 
efficiency.38  The  ratio  of  thrust  efficiency  to  energy  efficiency  was  asserted  to  describe  the  velocity  dispersion  in 
magnitude  and  direction.  This  ratio  will  be  derived  from  first  principles  in  the  following  section,  where  it  will 
become  apparent  that  the  ratio  is  equivalent  to  propellant  utilization  efficiency,  which  captures  losses  resulting  from 
plume  divergence,  incomplete  ionization,  and  the  production  of  multiply  charged  ions. 

The  term  propellant  utilization  efficiency  has  been  used  in  at  least  three  different  ways  in  the  literature. 

Most  commonly  as  (1)  the  ionization  fraction,  but  also  as  (2)  the  fraction  of  momentum  carried  by  ions18,  or  as  (3) 
the  ratio  of  output  moles  of  charge  to  the  input  moles  of  propellant39,  which  is  equivalent  to  the  product  of  the 
ionization  fraction  and  average  ion  charge.  While  the  definitions  of  voltage  and  current  utilization  are  more 
standardized,  a  consistent  description  of  propellant  utilization  has  not  emerged.  It  is  apparent  that  a  standard 
analytical  approach  is  needed  to  reconcile  the  terminology. 

The  analysis  of  voltage  utilization  and  current  utilization  have  been  adopted  in  several  performance  models, 
but  reduced  acceleration  potential  and  recycled  electron  current  have  not  been  associated  with  energy  efficiency. 
Recycled  electrons,  also  known  as  backflow  electrons  or  anode  electron  current,  were  characterized  as  a  dominant 
loss  mechanism  in  the  1960s.12  Hruby  et  al.  quantified  electron  recycle  fraction  as  a  loss  in  current  utilization 
efficiency,  and  accounted  for  loss  of  acceleration  potential  in  terms  of  reduced  voltage  utilization  efficiency. 5  A 
methodical  efficiency  analysis  by  Hofer  and  Gallimore  added  terms  to  account  for  the  utilization  of  mass  and 
charge.2  Several  other  studies  adopted  this  approach  and  expanded  the  analysis  to  include  factors  accounting  for 
plume  divergence  and  dispersion  of  the  jet  VDF.3'40' 41  These  and  other  recent  formulations  of  utilization  efficiencies 
do  not  establish  their  connection  to  fundamental  conservation  principles  and  Newton’s  second  law,  although  in  some 
instances  the  product  of  terms  is  remarkably  close  to  the  product  of  the  three  utilizations  set  forth  in  this  paper. 
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Decoupling  Energy  Efficiency  and  Propellant  Utilization 


Using  Newton’s  Second  Law,  steady  state  thrust  is  defined  in  Eq.  (1)  as  the  total  axial  momentum  along 
thruster  centerline.  Thrust  may  be  factored  into  the  product  of  total  mass  flow,  mass  weighted  average  velocity,  and 
momentum  weighted  average  divergence  using  the  definitions  shown  in  Eq.  (2),  (3),  (4),  and  (5). 
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The  term  v  ,  as  defined  in  Eq.  (3),  is  the  integral  over  all  velocity  space  of  the  VDF  at  azimuthal  position 
0.  All  quantities  are  assumed  steady  state  and  evaluated  at  constant  radius  from  the  exit  plane,  such  that  the 
quantities  to  be  integrated  have  units  of  sr'1. 

Anode  efficiency  is  expressed  in  terms  of  thruster  telemetry  and  measured  thrust  as  shown  in  Eq.  (6). 
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Further  analysis  of  thrust  efficiency  reveals  its  relationship  to  energy  efficiency  and  provides  insight  into 
loss  mechanisms  in  Hall  thruster  performance.  In  Eq.  (7)  and  (8),  thrust  efficiency  is  separated  into  the  product  of 
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energy  efficiency  and  propellant  utilization  using  the  definition  of  jet  kinetic  energy  (Pjet  =  Vi  m  <v2>„,),  Newton  s 
Second  Law,  and  the  mass  weighted  average  squared  velocity  as  defined  in  Eq.  (10). 
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Energy  efficiency  as  defined  in  Eq.  (11)  characterizes  the  conversion  of  input  anode  electrical  energy  to  jet 
kinetic  energy,  and  contains  all  information  about  losses  that  ultimately  appear  as  Joule  heating  of  the  thruster, 
radiation  from  the  jet,  and  frozen  ionization  losses  in  the  plume. 
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energy 
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The  quantity  <I>  as  defined  in  Eq.  (12)  is  a  mathematical  relationship  between  the  axial  momentum  that 
produces  thrust  and  the  jet  kinetic  energy.  It  contains  all  loss  information  associated  with  the  jet  vector  properties 
and  is  unity  for  100%  ionization  to  a  single  ion  species  where  all  thrust  vectors  are  directed  along  the  same  axis. 
This  ratio,  termed  propellant  utilization  efficiency,  accounts  for  losses  due  to  incomplete  ionization,  plume 
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divergence  and  non-uniformity  of  the  VDF  that  occurs  when  the  jet  is  composed  of  multiple  ion  species  with  widely 


varying  velocities. 


'  \ 2 

<  v  >„,<cos(e)>mv  _  2 

<]>  =  =  <cos(0)>ml,=<l>VDF<[>DIV  (12) 

<v2>ffl  <v2>m 

For  simplicity,  the  propellant  utilization  may  be  factored  into  two  parts.  The  term  <1>vdf  is  defined  in  Eq. 
(13)  and  quantifies  losses  from  non-uniformity  of  the  VDF  that  occur  when  the  jet  is  composed  of  multiple  particle 
species  with  widely  differing  velocities.  In  Eq.  (14),  the  term  <t>DIV  quantifies  jet  momentum  losses  from  plume 
divergence. 


<  v  >: 


<  v2  >„ 


(13) 


®DIV  — <COs(0)>"v 


(14) 


This  formulation  is  similar  to  the  focusing  efficiency  from  the  Kharkov  Aviation  Institute  methodology28, 
but  includes  the  effects  of  unionized  propellant  and  differentiates  between  mass  weighted  average  quantities  and 
momentum  weighted  average  quantities. 


Evaluating  Velocity  Distribution  Losses  in  Propellant  Utilization 


The  influence  of  non-uniform  velocity  distribution  on  propellant  utilization  is  expressed  in  Eq.  (15)  using 
the  definitions  in  Eq.  (4)  and  (10).  Knowledge  of  fly)  throughout  the  plume  is  required  for  a  rigorous  calculation  of 
this  influence  in  4>Vdf- 
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(15) 


To  illustrate  the  effects  of  non-uniform  VDF,  a  simple  trimodal  ion  population  with  neutral  xenon  is 
studied.  To  further  simplify,  the  velocity  of  the  neutrals  and  each  ion  species  is  approximated  with  a  delta  function 
distribution  of  velocities.  This  approximation  will  not  significantly  affect  conclusions  that  would  result  from  a  more 
complicated  model,  since  each  ion  species  is  created  within  a  narrow  zone  in  the  channel  and  the  plume  is 
predominantly  composed  of  Xe+1  ions.  The  value  of  4>Vdf  is  shown  in  Eq.  (16)  for  a  trimodal  ion  distribution  and 
evaluated  using  reduced  number  fractions  and  velocity  ratios  for  compact  notation.  Reduced  number  fractions  are 
defined  according  to  Eq.  (17)  and  (18)  where  fH=f,/fi,  f  1=fi/fi,  f\=f2/fi,  and  f  :=f  j/fr  Ratios  of  velocities  of  Xe°, 
Xe+1,  Xe+2,  and  Xe+3  to  the  velocity  magnitude  of  Xe+I  are  defined  by  y0=V(/lvil,  yi=Vi/lvjl,  y2=v2/lvil,  y3=V3/lvil. 


<  v2  >„ 


(v0f()+  Vifi+  V2f2+  y3f3)2  [  (yQ(i-fj)+  yifi  +  y2fa+  y3f3*)2 1 ..  f 

v0  f0+  V1  fl+  v2  f2+  v3  f3  L^O  (l-fi)+  yi  fl  +  y2  f2+  y3  f3  J 


where,  by  definition 

f0+(fl  +  f2  +  f3)=fo+fi=l 


(17) 


fi  +f2+f3  =1 


The  term  inside  the  brackets  of  Eq.  (16)  is  denoted  «.  The  factor  a  approaches  unity  as  f,  approaches  100% 
and  the  jet  composition  converges  to  a  single  ion  species.  In  a  rigorous  calculation  using  the  VDF  throughout  the 
plume,  «  may  be  evaluated  using  Eq.  (19). 


71/2  ' 

k  J  m(0)sin(0)  v  (0)d0 

^  VDF  _  1  V  -n/2 _ J 

f.  mf.  f  n!2 
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(19) 


The  factor  «,has  a  weak  dependence  on  ionization  fraction  for  f,  >0.70,  which  is  considered  a  lower  limit  on 
ionization  fraction  for  the  jet-mode  Hall  thruster  operating  regime.  Experimental  ion  species  compositions  from 
numerous  Hall  thruster  investigations  were  analyzed  using  the  assumption  of  a  trimodal  ion  population  with  delta 
function  velocity  distributions.  For  100%  ionization,  a  ranged  from  approximately  0.94-0.99,  and  is  shifted  upward 
2-3%  for  fj=0.70.46  The  trimodal  ion  composition  analysis  established  that  «  has  an  empirical  dependence  on 
average  charge  Q  for  delta  function  velocity  distributions,  which  is  expressed  in  Eq.  (20)  for  calculation  of  «  with 
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unknown  ionization  fraction.  Uncertainty  in  the  empirical  linear  relationship  is  conservatively  estimated  at  ±2%. 


Further  examination  of  «  in  Appendix  A  evaluates  variations  due  to  ionization  fraction  and  ion  species  composition. 


asl.14-0.13Q 


(20) 


where,  average  charge  Q  is  defined  by 

Q  =  f1*+2f2*+3f3* 


(21) 


Evaluating  Plume  Divergence  Losses  in  Propellant  Utilization 


The  momentum  weighted  average  plume  divergence  is  defined  in  Eq.  (22)  as  the  ratio  of  axial  momentum 
that  produces  thrust  to  the  total  momentum  exhausted  in  the  plume.  Momentum  losses  associated  with  plume 
divergence  may  be  calculated  with  knowledge  of  the  input  mass  flow,  measured  thrust,  and  the  mass  weighted 
average  velocity. 


®DIV  <  cos  ((()  >  mv 
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The  value  <cos(0)>"  has  been  used  in  previous  analyses  to  describe  plume  focusing"  ’  ’  ,  but  there  has  not 
been  a  consistent  method  to  calculate  the  effect  of  plume  divergence  on  thrust.  This  is  primarily  due  to  the  difficulty 
of  measuring  particle  velocity  throughout  the  plume.  Charge  divergence  in  the  plume  is  indicative  of  off-axis 
velocity  losses  in  thrust,  and  is  a  useful  alternative  for  experimental  characterization  of  performance  losses  due  to 
plume  divergence.  The  momentum  weighted  average  divergence  is  approximated  as  the  charge  weighted  average 
divergence  in  Eq.  (23),  which  enables  calculation  of  off-axis  cosine  losses  in  Eq.  (24)  using  the  ratio  of  axial  beam 
current  to  total  beam  current  as  measured  by  a  Faraday  probe. 
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(24) 
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Off-axis  cosine  losses  integrated  in  the  numerator  of  Eq.  (24)  quantify  the  axial  beam  current  that  generates 
thrust.  This  formulation  creates  a  method  where  the  plume  vector  loss  is  evaluated  in  a  scalar  form.  The  calculation 
of  <1>div  is  complicated  by  the  presence  of  charge  exchange  in  the  plume,  which  causes  artificial  increases  in 
measured  beam  current  at  high  angles  and  artificially  increases  divergence  losses. 

The  ratio  hi  (0)v(0)/j(0)  introduced  in  Eq.  (23)  reveals  the  difference  between  momentum  divergence  and 
charge  divergence.  This  ratio  is  evaluated  in  Appendix  B  for  variations  in  ionization  fraction  and  ion  species 
population.  In  Eq.  (25),  the  ratio  is  evaluated  at  off-axis  angle  0,  and  reduced  to  a  function  of  Q,  AV,  f;,  y0,  and  a. 


j(9)  l  F 


1/2 


\  1/2 


1/2 


=  |2Avf|  rj 


-,-1/2 


(25) 


The  term  inside  the  brackets  of  Eq.  (25)  is  denoted  y,  and  characterizes  differences  between  momentum 
weighted  divergence  and  charge  weighted  divergence.  The  factor  y  is  proportional  to  Q"1/2,  where  spatial  variation 
of  average  charge  in  the  plume  is  the  primary  source  of  discrepancy  between  momentum  weighted  divergence  and 
charge  weighted  divergence.  Significant  variations  in  average  charge  within  the  azimuthal  region  of  highest  beam 
current  may  alter  plume  divergence  calculations  by  less  than  3%.  Secondary  effects  evident  in  Eq.  (25)  include 
ionization  fraction  and  the  role  of  average  acceleration  potential  in  momentum.  Appendix  B  details  the  analysis 
used  in  Eq.  (25)  and  illustrates  variations  due  to  ionization  fraction  and  ion  population. 


Decoupling  Voltage  and  Current  Utilization 

Energy  efficiency  is  the  ratio  of  jet  kinetic  energy  to  anode  electrical  energy,  and  is  factored  in  Eq.  (26)  as 
the  product  of  voltage  utilization  and  current  utilization.  It  is  convenient  to  introduce  %,  which  is  defined  in  Eq.  (27) 
as  the  product  of  particle  ionization  fraction,  t't,  and  average  ion  charge,  Q,  for  a  plasma  with  three  ion  species. 
Although  x  cancels  in  the  energy  efficiency  term,  it  is  shown  in  Equations  (30)  and  (31)  to  reveal  the  physical  nature 
of  the  voltage  and  current  utilization  inherent  in  the  thrust  efficiency. 
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where,  by  definition 

X=  fj  +  2f2  +3f3  =  fi(f1*+2f24+3f3*)  =  f1Q  (27) 

The  ratio  of  average  particle  specific  kinetic  energy  to  average  ion  specific  kinetic  energy  in  Eq.  (26), 
<\2>J<\iom2>,  is  written  in  terms  of  f  and  reduced  species  mass  fractions  in  Eq.  (28). 
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The  bracketed  quantity  in  Eq.  (28)  is  very  close  to  unity  under  all  reasonable  conditions  of  Hall  thruster 
operation  due  to  the  low  value  of  y0  in  the  numerator,  i.e.,  y02(l-fi)=!lxlO"4  for  experimental  estimates  of  v, -0.03. 12 
Thus,  the  approximation,  <v2>m/<V;0,!S2>m  =  f,  is  accurate  to  better  than  one  in  ten  thousand. 

Equation  (29)  transforms  the  first  term  in  brackets  of  Eq.  (26)  into  an  expression  that  contains  explicit 
utilization  of  acceleration  potential  for  each  ion  species,  §V|/Va,  SV2/Va,  8V3/Va.  An  ExB  probe  measurement  of 
reduced  ion  mass  fractions  and  estimation  of  ion  acceleration  potentials  would  enable  a  voltage  utilization  to  be 
calculated  for  each  ion  species,  since  the  particle  kinetic  energy  is  proportional  to  charge  and  acceleration  potential 
(Vi  vz2=  z5Vz  &/J(). 
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(29) 


To  simplify  the  analysis,  all  ions  are  considered  to  be  created  in  the  same  zone  whose  length  is  small 
compared  to  the  acceleration  length,  such  that  AV=8Vi=8V2=8V3.  Previous  investigations  by  Kim43  and  King44 
found  that  species  dependent  energy  to  charge  ratios  varied  by  tens  of  volts  using  different  types  of  energy 
analyzers.  This  variation  energy  to  charge  ratio  is  less  than  10%  of  the  discharge  voltage.  According  to  Hofer45,  the 
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approximation  will  have  a  negligible  effect  on  accuracy  since  the  plume  is  predominantly  composed  of  singly 
ionized  xenon.  Thus,  the  most  probable  ion  acceleration  potential,  AV,  as  measured  with  a  retarding  potential 
analyzer  (RPA)  enables  the  first  term  in  brackets  of  Eq.  (26)  to  be  expressed  in  terms  of  the  average  voltage 
utilization  efficiency  and  %  in  Eq.  (30). 


-<v2 

<v->„,  2  w’ 

<v2  > 
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(f;+2f2*  +3f;)=f(i-p)Q  =  (i-p)X 


The  average  voltage  utilization,  (1-P)=  AV/Va,  compares  the  most  probable  potential  gained  by  ions  with 
the  applied  anode  potential,  and  is  unity  for  ionization  at  the  anode  face  where  ions  are  accelerated  through  the 
entire  anode  potential. 

The  second  bracketed  term  in  Eq.  (26)  is  transformed  in  Eq.  (31)  into  an  expression  containing  the  current 
utilization  efficiency  and  y_. 

m  F  ^  m(F/M)  f,Q  =  m(F/M)  x  1  =  ^Beam  1  =  (l~r)  (31) 

Ia  M  Ia  fjQ  \  1  \  1  1 


Current  utilization  efficiency,  (l-r)=  Ieean/Ia,  is  the  fraction  of  cathode  electron  flow  that  electrically  neutralizes  the 
accelerated  positive  ions  in  the  plume,  and  is  defined  as  the  ratio  of  ion  beam  current  to  discharge  current.  The  total 
ion  beam  current  exiting  the  thruster  is  given  by  lBeam=  X^1 

Electrons  recycled  to  the  anode  ionize  neutral  propellant  and  sustain  dissipative  plasma  processes.  The 
minimum  power  required  to  produce  a  given  trimodal  ion  distribution  is  mf;  (f:  i  f  +  e2f2  +  r.if<  )  such  that  the 
minimum  recycled  electron  fraction  required  to  sustain  ionization  may  be  expressed  by  Eq.  (32). 


(32) 


Power  from  recycled  electrons  that  is  not  used  for  ionization  is  lost  by  Joule  heating  and  amounts  to  (r-rmin)IaVa. 

Combining  the  derivation  of  energy  efficiency  with  propellant  utilization  efficiency  in  Eq.  (33)  yields  two 
useful  groups  of  experimental  parameters,  denoted  E,  and  E2,  which  may  be  written  using  the  preceding  utilization 
efficiencies.  The  parameters  are  formed  such  that  the  product  of  Ei  and  E2  is  equal  to  thruster  efficiency.  Equations 
(34)  and  (35)  show  how  these  quantities  are  calculated  based  solely  on  telemetry  and  thrust  measurements,  and 
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provide  insight  about  the  relative  magnitudes  of  the  individual  utilization  efficiencies  in  the  absence  of  plume 
measurements. 
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For  fixed  thruster  telemetry  inputs  Va  and  m  ,  the  dimensionless  experimental  parameters  E[  and  E2  separate 
changes  in  anode  efficiency  to  processes  related  to  variations  in  discharge  current  or  thrust.  The  quantity  E|~T2,  and 
is  a  function  of  propellant  utilization,  voltage  utilization,  and  X-  Experimental  parameter  Ei  relates  the  applied 
acceleration  potential  to  the  jet  vector  properties.  The  quantity  E2~Ia  \  and  is  a  function  of  current  utilization  and 
the  X-  Experimental  parameter  E2  relates  the  input  mass  flow  to  the  flow  of  charge.  The  inverse  of  E2  was 
developed  in  the  Russian  manual  on  stationary  plasma  engines  in  1989,  where  it  was  defined  as  the  exchange 
parameter.28  The  naming  convention  used  here  describes  the  exchange  of  applied  input  parameters  to  operational 
thruster  properties.  Thus,  E,  is  appropriately  termed  the  Voltage  Exchange  Parameter  and  E2  is  termed  the  Mass 
Exchange  Parameter.  While  ionization  and  acceleration  processes  are  closely  coupled,  the  form  of  the  experimental 
parameters  indicates  that  propellant  utilization  and  voltage  utilization  are  principal  in  the  formation  of  thrust.  In  the 
absence  of  diagnostics  for  the  determination  of  plasma  properties,  these  experimental  parameter  groups  allow  limits 
to  be  placed  on  acceptable  values  for  the  average  charge  Q,  ionization  fraction  f„  and  momentum  divergence  loss 

(D  46’47 

'A' rav¬ 
in  the  absence  of  thrust  measurements,  such  as  on-orbit  applications,  thrust  may  be  estimated  using 
telemetered  beam  voltage  and  beam  current  as  shown  in  Eq.  (36)  and  (37).  This  formulation  is  similar  to  the  one 
dimensional  analysis  used  to  estimate  thrust  of  the  SERT-II  ion  engines,  which  is  equivalent  to  Eq.  (37)  when  %  and 
<t>  are  unity.48,49,50'51'52'53  An  analogous  expression  was  used  to  estimate  thrust  of  the  NSTAR  ion  engines  on-board 
NASA’s  Deep  Space  1  by  including  factors  to  account  for  multiply  charged  ions  and  beam  cosine  losses.54  These 
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factors  are  based  on  the  formulation  of  thrust  by  Masek  et  al.  and  are  related  to  the  ratio  ( <!>/%) 1/2  as  shown  in  Eq. 


(37). 
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Experimental  Demonstration  of  Analytical  Methodology 

The  efficiency  analysis  is  demonstrated  using  experimental  results  from  the  NASA-173Mv2  laboratory 
model  Hall  thruster  operating  at  700  V  and  10  mg/s.2'45'55  The  investigation  consisted  of  thrust  measurements  and 
examination  of  the  ion  voltage  distribution,  plasma  potential,  ion  species  composition,  and  ion  current  density  in  the 
far-field  plume.  Experimental  results  are  listed  in  Table  1. 

Calculated  plasma  properties  in  Table  2  show  a  high  ionization  fraction  of  -98%  and  a  large  ion  beam 
current.  The  backflow  electron  current  not  applied  to  beam  neutralization  is  r=0.19,  whereas  the  minimum  recycled 
electron  fraction  required  for  propellant  ionization  is  rmin=0.02  of  the  total  discharge  current.  Losses  due  to  Joule 
heating  are  estimated  at  Ir-r^J^Va®!  100  W,  whereas  the  total  energy  losses  are  (l-qenergy)IaVaS!1500  W  of  the  total 
6700  W  discharge  power. 

Examination  of  the  propellant  utilization  shows  that  plume  divergence  is  a  dominant  performance  loss 
mechanism  for  this  operating  condition.  For  <J>DIV=0.88,  the  charge  weighted  plume  divergence  is  20°,  whereas  the 
commonly  reported  95%  charge  divergence  half-angle  is  significantly  higher  at  39°.  Losses  due  to  non-uniformity 
of  the  VDF  were  minimal  for  ,i>vni=0.97,  where  the  magnitude  of  a  was  estimated  using  the  experimental  ion 
species  fractions  with  delta  function  velocity  distributions. 

Thruster  utilization  efficiencies  and  anode  efficiency  are  shown  in  Table  3.  Anode  efficiency  is  calculated 
using  two  methods,  (1)  using  measured  thrust  and  telemetry  and  (2)  using  plume  diagnostics.  Both  methods  resulted 
in  anode  efficiency  of  66%.  The  overall  uncertainty  in  thruster  efficiency  based  on  thrust  and  telemetry  data  is  2%, 
whereas  the  uncertainty  in  efficiency  using  plasma  diagnostics  is  -15%. 
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The  uncertainty  in  thrust  and  plume  measurements  were  propagated  into  calculations  of  plasma  properties 
and  utilization  efficiencies.  Experimental  errors  associated  with  Faraday  probe  scans  manifest  in  both  current 
utilization  and  propellant  utilization,  and  translate  into  a  large  degree  of  uncertainty  when  calculating  thruster 
efficiency  using  plume  diagnostics.  As  distance  from  the  exit  plane  increases,  charge  exchange  (CEX)  collisions  in 
the  plume  lead  to  a  decrease  in  core  ion  beam  current  and  an  increase  in  computed  beam  divergence.  These  effects 
artificially  decrease  efficiency  as  calculated  with  plume  diagnostics  and  increase  measurement  uncertainty.  Near¬ 
field  measurements  close  to  the  exit  plane  would  minimize  the  effects  of  CEX  collisions. 

The  far-field  Faraday  probe  measurements  on  the  NASA-173M2  were  taken  at  98.5  cm  downstream  of  the 
exit  plane,  leading  to  estimated  uncertainty  in  axial  and  total  beam  current  of  -10%.  If  the  uncertainty  in  calculated 
beam  current  calculations  were  reduced  to  5%,  the  overall  uncertainty  in  efficiency  calculated  with  plume 
measurements  would  decrease  from  15%  to  8%.  This  highlights  the  importance  of  high  fidelity  measurements  with 
minimal  uncertainty  for  calculating  plasma  properties  using  the  analytical  technique. 


Table  1:  Thrust,  Telemetry,  and  Far-field  Plume  Measurements 


T  [mN] 

298  ±  1% 

Va  [V] 

700.3  ±  0.05% 

la  [A] 

9.62  ±  0.2% 

m  [mg/s] 

10.0  ±  1% 

f\ 

0.92  ±  0.04 

0.07  ±  0.02 

f*3 

0.01  ±0.01 

AV  [V] 

675  ±  10 

I-Beam  [A] 

7.8  ±  10% 

1  Axial  [A] 

7.3  ±  10% 

Table  2:  Calculated  Plasma  Properties 

X  =  Ib  /  ( A  &/JC) 

1.07  ±  10% 

Q  =  f\  +  2f*2  +  3f*3 

1.09  ±3% 

f,=  x/Q 

0.98  ±  10% 

«  =  (0.03(1  -fd+f*  !+V2f*  2+V3f*  3)2/(0.00 1  ( 1  -f;)+Q) 

0.99  ±  2% 

®VDF  =  «  f 

0.97  ±11% 

®DIV  =  (lAxial  /  teamf 

0.88  ±  10% 

X  =  cos'1  V  ®MV 

20  °±  10% 

95%  Charge  Divergence  Half-Angle 

39° ±  10% 

rmm  =  m  fi  (e  / !+£/ 2±£3f* 3)/(Va  Ia) 

0.02  ±11% 

16 


Table  3:  Calculated  Thrust  Efficiency  and  Utilization  Efficiencies 


0.87  ±  2% 

E2=(  m  /  Ia)  (  &/JC) 

0.76  ±  1% 

(1-(3)  =  AV/Va 

0.96  ±1% 

( 1  -l)  =  I-Beam  /  I 

0.81  ±  10% 

henergy  —  (1"P)  (1+) 

0.78  ±  10% 

®  =  ®VDF  ®DIV 

0.85  ±  15% 

Hthrust  =  (1-p)  (1-r)  <E>  (Analytical) 

0.66  ±  15% 

Thrust  =  ‘/2  T2/  (  rh  Va  Ia)  (Telemetry) 

0.66  ±  2% 

Summary 

A  standardized  Hall  thruster  performance  architecture  was  developed  that  factored  anode  efficiency  as  the 
product  of  terms  representing  the  utilization  of  thruster  input  parameters,  including:  (1)  voltage  utilization,  (2) 
current  utilization,  and  (3)  propellant  utilization  efficiency.  The  analysis  decomposed  anode  efficiency  from  first 
principles,  and  formulated  the  relationship  such  that  losses  associated  with  energy  conversion  are  analyzed 
separately  from  losses  associated  with  jet  vector  properties.  Energy  efficiency  encompassed  performance  losses 
resulting  from  Joule  heating  and  ionization  processes,  and  is  defined  as  the  product  of  current  utilization  and  voltage 
utilization.  Current  utilization  can  never  be  unity  due  to  a  finite  flow  of  electrons  that  are  recycled  to  the  anode  to 
sustain  ionization  processes.  Propellant  utilization  efficiency  incorporates  losses  from  plume  divergence  and  the  ion 
jet  composition,  and  is  unity  for  100%  ionization  to  a  single  ion  species  whose  velocity  vectors  are  directed  along 
the  same  thrust  axis.  The  momentum  weighted  divergence  loss  term  (<cos(0)>2mv)  in  propellant  utilization  contains 
all  jet  vector  losses.  The  proposed  analysis  is  compared  to  past  performance  methodologies,  with  similarities  to 
previous  architectures  formulated  by  Masek  et  al.  in  the  ion  engine  community  and  Belan,  Kim,  Oranskiy,  and 
Tikhononv  from  the  Kharkov  Aviation  Institute.21'28 

A  case  study  was  presented  with  experimental  performance  and  far-field  plume  measurements  for  a 
laboratory  model  Hall  thruster.  Thrust  measurements  resulted  in  anode  efficiency  of  66+2%,  compared  to  66+25% 
when  efficiency  was  calculated  from  plume  measurements  using  the  analytical  model.  The  performance 
methodology  reveals  low  current  utilization  and  high  plume  divergence  as  the  primary  performance  drivers  for  this 
operating  condition.  Although  uncertainty  associated  with  Faraday  probe  current  density  measurements  in  the  jet 
were  significant,  the  numerical  agreement  between  the  two  techniques  supports  the  proposed  methodology.  Fidelity 
of  the  performance  architecture  and  the  nature  of  analytical  factorization  permit  the  extraction  of  additional  plasma 
parameters  and  the  placement  of  bounding  limits  on  discharge  properties. 
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APPENDIX  A:  The  VDF  Factor  «Tor  a  Trimodal  Ion  Distribution 


The  <t>VDF  component  of  propellant  utilization  is  often  absent  in  textbook  equations  where  scalar  energy 
quantities  are  compared  to  squared  vector  momentum  quantities,  since  <t>VDF  >  0.99  for  shifted  Maxwellian 
distributions  common  to  chemical  propulsion  jets.  To  illustrate  the  dependence  of  <J)  on  the  properties  of  a  jet  with 
non-uniform  velocity  distribution,  a  simple  trimodal  ion  distribution  with  neutral  xenon  is  studied.  Consider  a  jet 
that  contains  Xe°,  Xe+1,  Xe+2,  and  Xe+3  with  mass  fractions  of  f0,  f),  f2,  and  f3,  respectively,  such  that  f0+fi+f2+f3=l. 
To  further  simplify,  the  velocity  of  neutrals  and  each  ion  species  is  approximated  with  a  delta  function  distribution 
of  velocities,  specified  by  v0,  Vi,  v2,  and  v3.  This  approximation  will  not  significantly  affect  conclusions  that  would 
result  from  a  more  complicated  model,  since  each  ion  species  is  created  within  a  narrow  zone  in  the  channel  and  the 
plume  is  predominantly  composed  of  Xe+1  ions.  The  term  <t>VDF  is  shown  in  Eq.  (A-l)  for  a  trimodal  ion 
distribution,  and  the  extension  to  include  higher  charge  states  is  a  simple  exercise. 


<v>»,  _  (v0f0+  vlfl+  v2f2+  v3f3)2 

<  V"  >m  V02f0+  V,2f,+  V22f2+  V32f3 


(A-l) 


For  convenience  and  compact  notation,  the  ion  number  fraction  is  defined  by  fi=f1+f2+f3,  where  f0+f;=l. 
Reduced  number  fractions  are  defined  by  f  n=f{/fi,  f  1=f1/fi,  fT2=f2/f;,  and  f  3=f3/f;,  such  that  f  1+f*2+f*3=l.  Ratios  of 
velocities  of  Xe°,  Xe+1,  Xe+2,  and  Xe+3  to  the  velocity  magnitude  of  Xe+I  are  defined  by  y0=Vo/lvil,  v i = v ,/l v [ I, 
y2=v2/lvil,  y3=v3/lvil.  This  enables  0Vi)[  to  be  written  in  Eq.  (A-2)  as  a  function  explicit  in  I’  and  the  reduced  ion 
mass  fractions  t)  .  f2*,  and  f3  . 


=f.  (yod-fjH  yi fi*+  y2f2+  y3f.1l2  =f. . 
1_y02(1-fi>|-  y  1 2 f T  +  y22f2+  y32f3 .  1 


The  term  inside  the  brackets  of  Eq.  (A-2)  is  denoted  «  for  a  trimodal  ion  composition.  The  factor  » 
approaches  unity  as  f,  approaches  100%  and  the  jet  composition  converges  to  a  single  ion  species.  Figure  A-l 
shows  the  relationship  between  «  and  average  ion  charge,  Q,  for  a  trimodal  ion  population  with  100%  ionization. 
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Reduced  ion  mass  fractions  f  b  f  2,  and  f  3  are  determined  from  several  experimental  studies  of  the  three  ion  species 
fractions,  depicting  a  range  of  thrusters,  facilities,  and  operating  conditions.45’56’57'58'59'60  To  illustrate  the  variation  in 
»  using  Eq.  (A-2),  ion  velocity  ratios  of  each  species  are  calculated  for  the  idealized  case  of  ion  creation  in  the  same 
zone  whose  length  is  small  compared  to  the  acceleration  length.  This  approximation  is  consistent  with 
measurements  by  Kim43  and  King44,  who  found  that  species  dependent  energy  to  charge  ratios  varied  by  tens  of 
volts  using  different  types  of  energy  analyzers.  As  pointed  out  by  Hofer45,  this  assumption  will  have  a  negligible 
effect  on  the  ion  velocity  ratio  estimates  since  the  plume  is  predominantly  composed  of  singly  ionized  xenon.  Thus, 
the  ion  species  kinetic  energies  are  proportional  to  their  charge,  such  that  the  velocity  ratio  magnitudes  are  lyil=l, 
ly2l=V2,  and  ly3l=V3.  The  velocity  of  neutrals  is  at  least  30  times  smaller  than  that  of  Xe+1,  resulting  in  lyol^O.OS.42 


Average  Charge,  Q 

Fig.  A-l  Dependence  of  »on  plume  ion  composition.  A  complete  set  of  constant  f( ,  IV  ,  f3*  lines  are  shown  for  the  f  =  1.0 
manifold.  The  factor  «ds  calculated  from  ion  species  fractions  reported  in  the  Hall  thruster  literature  for  delta 
function  velocity  distributions,  including:  NASA-173Mv2  (2.5-9  kW)45,  ATON  (~1  kW)57,  BHT-600  (0.6  kW)57, 
SPT-100  (1.2  kW)58,  P5  (5  kW)59,  and  the  BHT-4000  (3  kW)60. 


The  factor  a  has  a  weak  dependence  on  ionization  fraction  for  f,  >0.70,  which  is  considered  a  lower  limit  on 
ionization  fraction  for  all  jet-mode  operating  regimes.  The  non-linear  sensitivity  of  «  to  ionization  fraction  is 
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illustrated  in  Figure  A-2  for  manifolds  of  £=0.2  and  £=  1 .0.  The  region  from  1.0<Q<1.5  is  shaded  for  0.7<  £<1.0, 


which  brackets  the  magnitude  of  «  from  0.93  to  1 .03  for  most  Hall  thruster  operation. 


Average  Charge,  Q 

Fig.  A-2  Variation  in  a  due  to  ionization  fraction,  average  charge,  and  ion  species  population  for  a  trimodal  ion 

composition  with  delta  function  velocity  distributions.  Ionization  fraction  manifolds  are  displayed  for  £  =  1.0 
and  £  =  0.2.  Dashed  lines  of  constant  If  =  0  and  f3*  =  0  are  shown  for  £  =  0.9  to  0.3  in  increments  of  0.1. 


The  shaded  region  in  Figure  A-2  is  magnified  in  Figure  A-3,  and  the  magnitude  of  »is  evaluated  for  0.7 < 
£<1.0  to  illustrate  the  range  of  «  for  the  survey  of  Hall  thruster  ion  species  compositions.  For  100%  ionization,  « 
ranged  from  approximately  0.94-0.99,  and  is  shifted  upward  2-3%  for  £=0.70.  The  studies  establish  an  empirical 
linear  relationship  between  «  on  Q,  which  is  expressed  in  Eq.  (A-3).  Uncertainty  in  «  using  the  empirical 
relationship  is  conservatively  estimated  at  +2%. 

a  =  1.14-0.13Q  (A-3) 
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Average  Charge,  Q 

Fig.  A-3  Range  of  a  from  the  survey  of  Hall  thruster  Ion  species  compositions  for  a  trimodal  ion  population  with  delta 

function  velocity  distributions.  The  f|  =  1.0  manifold  is  displayed  along  with  dashed  lines  of  constant  f3*  =  0  for  I- 
=  0.7,  0.8,  and  0.9.  The  magnitudes  of  a  are  plotted  for  fj  =  0.7,  0.8,  0.9,  and  1.0. 
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APPENDIX  B:  The  Divergence  Factor  y  for  a  Trimodal  Ion  Distribution 


The  momentum  weighted  average  divergence  is  approximated  as  the  charge  weighted  average  divergence 
in  Eq.  (B-l),  which  enables  calculation  of  off-axis  cosine  losses  using  the  ratio  of  axial  beam  current  to  total  beam 
current  as  measured  by  a  Faraday  probe. 


v(e)>L= 


m(e)v(en . 

T5T  ' 


m  z. 

j(0)cos(0)sin(0)d0  2tt  J  j(0)cos(0)sin(0)d0 


III  z. 

J  j(0)  sin(0)d0 


=<cos(0)>7  (B-l) 


Analysis  of  the  ratio  m  (0)v(0)/j(0)  introduced  in  Eq.  (B-l)  characterizes  the  properties  that  cause 
differences  between  momentum  weighted  divergence  and  charge  weighted  divergence.  Consider  the  trimodal  ion  jet 
from  Appendix  A  with  mass  fractions  of  f0,  fh  f2,  and  f3,  respectively,  such  that  t'n+fj+tj+l^  fo+fpl.  Reduced  mass 
fractions  are  defined  by  f  n=f0/fj,  f  1=f1/fi,  f  2=12/!,,  and  f  3=f3/fi,  where  f  i+f  2+f  3=1.  Each  particle  species  is 
approximated  with  a  delta  function  velocity  distribution  and  velocity  ratios  are  calculated  for  the  idealized  case  of 
ion  creation  at  the  same  location,  such  that  each  ion  is  accelerated  through  the  same  potential.  As  a  result,  the  ion 
species  kinetic  energies  are  proportional  to  their  charge  and  the  velocity  ratio  magnitudes  are  lyjl=l,  ly2l=V2,  and 
ly3l=V3.  The  ratio  111  (0)v(0)/j(0)  is  evaluated  at  off-axis  angle  0  and  formulated  in  terms  of  AV,  Q,  a,  f,,  and  y0  in  Eq. 
(B-2),  (B-3),  (B-4),  and  (B-5). 


iiv(e)v(e)  _  iiv(e)v(e) 
j(0)  m(0)tj(0)Q(0)^ 


fjmfeyo  +  jjVi  +  f2*y2  +  f3*y3) 


f,  +  2f,+3f. 


m(9)T(9)  _ 

M 

(fdyo+fiVi+f^yz+fsVa) 

2  A 

1 

1/2 

1 

(y02fd+  y,2fi*+  y22f2*+  y32f31 

j(e) 

V1 

T 

jo2fo+  yi2fi>+  y22f2+  y32; 

P” 

f3  J 

Q 

m(9)v(9)  ,  |  M  \,  Vn  (y02(l  —  fj )+  f*+  2f*+  3f3j  1/2 
j(9)  Q  J 


rir(9)T(9) 

j(e) 


(B-2) 

(B-3) 

(B-4) 

(B-5) 


The  term  inside  the  brackets  of  Eq.  (B-5)  is  denoted  y,  and  describes  differences  between  momentum 
weighted  divergence  and  charge  weighted  divergence  associated  with  average  charge  Q,  ionization  fraction,  and  ion 
species  population.  The  velocity  ratio  normalization  term,  lvjl=(2  AV  ■’W/Ji) 1  2,  shows  the  momentum  dependence  on 
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acceleration  potential.  The  quantity  y  is  further  reduced  in  Eq.  (B-6)  and  expressed  in  terms  of  <1>vdf  and  X-  F°r 
experimental  estimates  of  lyol®  0.03, 42  the  value  y02(  1  -fi)~ lxl 0"4.  Ionization  fraction  and  ion  species  population  are 
shown  to  have  a  negligible  effect  on  y  in  Figure  B-l.  Thus,  the  approximation  y  =  (a/Q)1/2  is  better  than  one  in  ten 
thousand. 


(B-6) 


Average  Charge,  Q 

Fig.  B-l  Variation  in  y  due  to  ionization  fraction,  average  charge,  and  ion  species  population.  Ionization  fraction 
manifolds  are  shown  for  f  =  1.0  and  f,  =  0  with  lines  of  constant  fi*=0,  IV  =0,  and  f3*  =  0,  and  f3*  =  0. 


Variations  in  the  magnitude  of  Q  alter  y,  but  the  azimuthal  location  of  variations  in  Q  is  the  dominant  factor 
causing  differences  between  momentum  weighted  divergence  and  charge  weighted  divergence.  This  effect  is 
illustrated  in  Figure  B-2,  which  shows  a  representative  distribution  of  current  density  and  beam  current  as  a  function 
of  azimuthal  position  for  a  laboratory  Hall  thruster  operating  at  6  kW.  The  location  of  peak  beam  current  in  the 
plume  is  approximately  8  degrees  from  thruster  centerline,  and  variations  in  Q  within  the  azimuthal  range  of  full 
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width  half  maximum  (FWHM)  will  have  the  greatest  effect  on  divergence.  For  the  beam  current  in  Figure  B-2,  this 
range  is  approximately  +3  to  +21  degrees.  Near  and  far-field  plume  measurements  of  the  SPT-100  and  BHT-200 
found  that  the  ion  species  fraction  did  not  significantly  change  within  +15  degrees  off  thruster  centerline,  and  the 
fraction  of  Xe+2  sharply  increased  by  approximately  5%  at  +20  degrees.  43'58'61  These  changes  in  ion  species 
fractions  have  a  negligible  effect  on  the  value  of  <cos(0)>mv,  thus  enabling  the  approximation  of  equivalence 
between  momentum  weighted  divergence  and  charge  weighted  plume  divergence  as  measured  with  a  Faraday  probe. 


-90  -75  -60  -45  -30  -15  0  15  30  45  60  75  90 

Azimuthal  Position  [Degrees] 


Fig.  B-2  Distribution  of  beam  current  and  current  density  at  1  meter  radius  as  a  function  of  azimuthal  position  in  the 
plume  of  a  laboratory  Hall  thruster  operating  at  6  kW. 
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